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Abstract 



We investigate a gauged matrix model in the large limit which is closely 
related to the superconductor fluctuation and the flux lattice melting in two 
dimensions. With the use of saddle point method the free energy is expanded 
up to eighth order for the coupling constant g. In the case that the coefficient 
of quadratic term of the Ginzburg-Landau matrix model is negative, a critical 
point (7 = is obtained in the large limit and the relation between this 
phase transition and the 2D flux lattice melting transition is discussed. 
PACS numbers: 74.40. +k, 74.60.-w 
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I. INTRODUCTION 



In the study of the fluctuation of superconductors and the melting transition of the 
Abrikosov flux lattice in a strong magnetic field, Ginzburg-Landau (GL) model is especially 
useful with the lowest Landau level approximation. This approximation is valid in the 
region very near to ifc2 and neglects all the contribution from upper Landau levels except 
the lowest one, and reduces the effective dimension of the system by two due to the Landau 
quantization perpendicular to the magnetic field axis. 

The perturbational studies |T[]- |^ for GL free energy up to high orders are unable to 
find a flux lattice melting transition in two dimensions. However, the comparison of the 
low temperature expansion to the high temperature expansion for the free energy gives the 
estimation of the melting transition point . On the other hand, the numerical studies such 
as Monte Carlo simulations show that there exists the first order transition in a finite tem- 
perature well below the mean field superconducting transition point Hc2 0- 0- Although 
both analytical and numerical calculations agree quite well in derivation of the statistical 
amounts, the nature of the melting transition of flux lattice is not yet fully understood |T0| - 



The first order transition has been suggested by the renormalization group analysis 



|T^, but two dimensions are far from the valid region around six dimensions. 

In this study we generalize the GL Hamiltonian to a gauged matrix model in which the 
order parameter is expressed by a complex N x N matrix ||13|. The matrix model in the 
large limit attracts recently much theoretical interest since it has a close relation to the 
string field theory and has been studied in many fields such as 2d quantum gravity 



coupled to matter field [jT5[, mesoscopic fluctuations [0 and electron correlations [0. In 



several matrix models the exact solutions have been obtained [T^ and these solutions give 
the clues to the perturbative analysis of other unsolved matrix models. It is known that 
the matrix model has a phase transition in the large A^ limit when the coefficient of the 
quadratic term is negative or in the lattice gauge theory [|^- This may be true for our 
gauged matrix model and we are interested in how this phase transition is related to the 
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superconducting flux lattice melting. Usual Ginzburg-Landau model corresponds to = 1 
case, but we generalize the order parameter to N x N complex matrix and take the large 
N limit. This large limit should be distinguished from the large case of the A^- vector 



model p2| , |23| , which becomes equivalent to the Hartree- Fock approximation and has no 
phase transition in two dimensions. 

In the usual GL model, the perturbation series becomes asymptotic expansions and one 
needs Borel summation for such divergent series. In the large A^ limit of GL matrix model, 
the perturbation series about g is convergent, and the precise analysis of the free energy 
becomes possible. Thus it is natural to consider the solution of the GL matrix model in the 
large A^ limit as a flrst approximation of the 2d superconductor phase transition. We flnd 
indeed a phase transition point gc/ot\ = 0.07 for the large A^ limit of the GL matrix model 
and this point corresponds to y = —2.7 in term of the reduced relative temperature. In the 
previous paper 0], we evaluated the flux lattice melting temperature for 2d superconductors 
(A^ = 1) as ?/ = — 10 by the analysis of the perturbation series of the usual GL free energy. 
Although the phase transition for A^ = oo occurs in considerably higher temperature than the 
usual GL model of A^ = 1, we consider this difference may result from the flrst approximation 
for 2d superconductors. The improvement of the evaluation of this phase transition point in 
the gauged matrix model is suggested by taking the higer order terms in 1/A^^ expansion. 

This paper is organized as follows : In section II, a gauged matrix model with the 
interaction of Ginzburg-Landau type -^Tr(M*M)^ is introduced in the lowest Landau level 
approximation and the free energy is obtained as the power series of coupling constant g 
using the saddle point method. In section III, the critical point of the free energy is obtained 
for the case the coefficient of TrM*M term is negative. This phase transition is investigated 
in nature and in section IV the relation to the melting transition of flux lattice is discussed. 
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II. GAUGED MATRIX MODEL IN TWO DIMENSIONS 



In the Ginzburg-Landau model, the order parameter is a complex field and can be 
expanded with the Landau levels. If the field is near the critical point Hc2, the lowest 
Landau level approximation, which neglects all contributions to the order parameter from 
the upper Landau levels and takes into account only the lowest one, is justified for the strong 
magnetic field due to the absence of mixing between Landau levels. We consider this strong 
magnetic field case and generalize this complex order parameter to a complex matrix 4>ij. 

The Hamiltonian is given by 

where is a rank N complex matrix and /i — x,y. We denote the charge of superconductor 
by e, and this charge e is of course twice of the charge of a single electron. Here the abehan 
gauge field ^4^ is a vector potential of a magnetic field, a and /3 are usual GL parameters. We 
choose the Landau gauge A = (0, Bx) . The complex matrix (p is written by the projection 
to the lowest Landau level as. 



/pB\ ^I'^ 



2 eB^ . 



(11.2) 



where Mjj is a. N x N complex matrix and we put m — 1 and h — 1. Then the Hamiltonian 
is rewritten as 

P feB\'/' 



X exp 



2eB 4 

X 5,,+,„,3+,,IV [M*{q,)M{qs)M*{q,)M{q,)] (II.3) 

where = a + ^ is related to the reduced temperature [T — Tc{B)]/Tc{B). The free 
energy F is obtained 

F = -^lnZ (II.4) 



dMe-^. (II.5) 



The perturbation series for the free energy is expanded by a new variable defined by 

This is the same as the expansion parameter in iV = 1 case 1^0] . Since we have generahzed 
the GL order parameter to a matrix variable, the perturbation series has a A^-dependence. 
In the diagrammatic language, the planar diagrams are obtained in the large limit. We 
have nonlocal interaction in Eg. ([II. 3D and have to perform Gaussian integration, which is 
equivalent to count the number of Euler path T of each diagrams [@]. Thus there appears 
additional factor 1/T to the combinatorial factor. In the large N limit, we have to select 
only planar diagrams among various terms which have been worked before by considering 
A^- dependence |l|]. 

Here we follow the new calculational method instead of the selecting planar diagrams as 



proposed in ref. [|T^. Except the factor T of the number of Euler path, the combinatorial 
factor of each diagram in the planar limit becomes the same as one matrix model. There is 
no difference between a complex matrix model and a Hermitian matrix model in the leading 
order except a trivial factor 2. It has been recognized that the renormalized expansion 
simplifies remarkably the diagramatic expansion for matrix models : one needs to consider 
only the irreducible diagrams to obtain the perturbation series of the free energy. This 
renormalized expansion can be applied successfully to our case, since the extra factor 1/T 
due to the Euler path is factorized. This factorization is easily understood for any diagrams 
by the definition of Euler path. Therefore, we follow the procedure of the renormalized 



expansion method [T^, which gives a remarkably simplified method for the perturbation 
expansion. 

We introduce the equivalent 2A^^-real vector model which is expressed by the 2A^^ di- 
mensional real vector field r and the effective Hamiltonian for large- limit is written as 

1 oo 

^efr = ^ + ]^EA/^'' (11.7) 



where 

x = r^ = (TtM*M) . (II.8) 

Note that TrM*M is the sum of the square of the absolute value for the matrix element. 
We integrate out the angular variables of this 2N'^- dimensional coordinates by keeping the 
radial part |r|. We choose the appropriate coefficients fk so as to make the free energy in 
the large limit becomes the same as that of the original gauged matrix model defined by 
Eq.( [11.3|) . This coefficient fk turns out to be determined by the irreducible diagrams in the 



original matrix model. 

By the saddle point method, the free energy for the Hamiltonian Eq. ( |II.3|) becomes in 
the large limit, 

fp °° 

-^ = x + J2fkg''x''-lnx (II.9) 

k=l 

where we have replaced x — ^ N^x. The saddle point equation is derived as 

= ^ + E mK)9'^'' -1 = 0. (11.10) 
This relation is also expressed simply by 

x = l-2gy (11.11) 

where we define 

The relevant irreducible diagrams in each order of g is reduced much in number and in 
Table I we give these numbers of diagrams. The self-energy part is completely renormalized 
into the quantity a; =< Tr M*M > and only irreducible diagrams should be taken into 
account. In other words, our irreducible diagrams do not contain the self-energy diagrams. 
We briefiy explain how we obtain the irreducible planar diagrams : we generate n-th order 
planar irreducible diagrams from {n — l)-th order ones by cutting two lines and add one 
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vertex. The choice of two hnes to be cut is restricted by the condition that generated next 
order diagrams should also be planar, and the directions of the lines are uniquely determined 
also by this condition. Firstly we make the generated graphs undirected and drop all the 
isomorphic graphs obtained by this process. Then for the remaining graphs, we choose 
properly directed graphs and calculate the combinatorial factor by counting the number of 
isomorphisms of the graphs. The number of Euler path of a given graph is determined by 
evaluating the determinant of the adjacency matrix 0. In Fig.l we show planar irreducible 
diagrams to eighth order, where the direction of the line (arrow) is not shown and the line 
is denoted simply by a sigle line instead of double lines. 

Up to the eighth order we obtain the perturbation series as 



where the coefficients are evaluated solely from the irreducible diagrams. 

It may be instructive to see how this renormalized expansion method gives an efficient 
result for the case of one matrix model where the exact solution is known : We now con- 
sider the case without a magnetic field in d = dimension. In the one matrix model, the 
Hamiltonian becomes with Hermitian matrix M 



from which we obtain the perturbation series for the free energy in the large N limit with 
the same irreducible diagrams. (Here we omit the factor 1/T for the Euler path.) Since M 




(11.13) 




(11.14) 



is Hermitian in Eq. (|II.14|) , there appears a difference of a factor 2 compared to the complex 



matrix case except the number T of Euler path. 



F \ 1 32 
— = — lna; + -a; + 2gx^ - 2o V H g^x^ - 96o V H 




(n.l5) 



The saddle point equation is given by 
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X ( dF\ 11 9 4 OR 



(11.16) 



From Eg. ( |li.l5| ) we obtain for y = d{F/N'^)/dg, 

y = 2x'^ - Agx^ + 32^V _ 334^3^8 ^ Q^^Ay ^^[.17) 

In the iteration for small g with Eg. ( [11. 17]) , is reexpressed as 

x^ = \y + \9v' - gV + ^qY (n.is) 



By the investigation of the ratio of the coefficients in the series Eg. (|II.18| ), we find that the 



ratio Rk = Ck/ck-i is given by Rk = — 12 + 30/A; {k > 3). Then we have a following closed 
eguation, 

1 



^ 108 



;i + Ugy f - 1] + f • ("-19) 



From this result, we can determine the coefficients of the free energy in Eg.(|LT5) up to any 
order of g for the one matrix model. 

We now come back to the case with a strong magnetic field. Including the factor 1/T 
of Euler path, we are able to check the result of Eg. (|II.13|) comparing with Eg. (|II.15|) . We 
repeat the same procedure as the one matrix model. From Eg. ([11.13| ) the saddle point 
eguation becomes 

I, (f^) = -1 + X + Agx^ - 2,V + - 1^,V + 0(/) 

= (11.20) 

We rewrite this eguation by the iteration for small g as 

9 1120 n 

x = l-Ag + Mg^ - ^g"" + • • • (11.21) 

Using the relation of Eg. ( [11. 11] ) (note y is defined hy y = d{F/N'^)/dg), the free energy is 
obtained in the large N limit, 
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Eiitim . 2, _ 11,. + 62H<,3 _ 58511 272 

^ 2 ^ 9^ 15^ 525^ 

-77001^/ + 993805 ^'^^J^^^ _ 1.351344723 x lO'g' (11.22) 

5005^ 33948915^ ^ ^ ^ 

We correct here the result of the previous calculation of order g''^ and in ref. [|T^, where 
small deviation was represented. It is remarkable that we easily obtain the series expansion 
by considering rather small diagrams. 

III. THE CRITICAL POINT 



Our main object is to apply the matrix model to the investigation of the flux lattice 
melting transition of a superconductor in a strong magnetic field for two dimensions. Here- 
after we consider the gauged matrix model with negative mass < and positive coupling 
case g > . 

From the study of the large limit of the two dimensional U (N) lattice gauge theory 
1^,20| and one matrix model |21|, it is known that there is a critical point of g and the 



free energy shows different behaviors in the small coupling {g < g^) and the strong coupling 
{g > gc) regions. Also the third derivative of the free energy with respect to the temperature 
becomes discontinuous at this point and third order phase transitions are observed in these 
models. This phase transition is equivalent to the freezing of the saddle point value of x. In 
the equivalent 2iV^- vector model representation, x is frozen below this critical point as 

X = ^ (III.1) 
9 

where s is a certain negative value to be determined. 

We again go back to the relation between x =< TyM*M > and y =< Tr(M*M)^ >. 



From Eq.( [11.13| ) y is expanded by 



y = = 2x' - gx' + Ig'x' - + • • • (III.2) 

This relation remains true even when we introduce the negative as- We reexpress this 
perturbation series as 
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2 -L -^2 2 '^ss 299 A A 8249 ^ c 

4624511 83706754319 . . 

H 9 V q y 

6589440^ ^ 49662412800^ ^ 

26766869658031714037 o « ,,,, , 

+ 5471161812032716800 ^ 

For the discussion of the region below the mean field transition temperature, as < 0, we 
write the free energy with an explicit a^-dependence in the term of order x, 

F = -Inx + obx + 2gx'^ -^g'^x^ -\ (111.4) 



In this case, we change the definition of g from the previous one Eq. (|II.6| ) to 



g = (III.5) 

The saddle point equation becomes, 

asx = 1 — 2gy (III. 6) 

Thus we have for the temperature as < 0, with Eq.( [lll.l] ) 

9 - sal /TTT 

yy = — ^ — ni.7 

As shown in ref. |jl^ , the critical point is obtained exactly for the one matrix model from 
the exact equation of Eq. ( [1I.19| ) by insertion of Eq. (|III.l|) and x = 1 — Agy. For our present 
model, the exact equation is unknown. Thus we approximate Eq.( |lII.3D , which corresponds 
to Eq.( [11.19| ) for one matrix model by the Fade form, and repeat the same analysis which 
determines the critical point. Approximated Fade form of Eq. ( [111.3]) becomes 

2 I + cigy + C2{gyy + ■ ■ ■ + Cg{gy)'i 
We put Eq.( |III.7|) into this equation then we have 

where z = —sa\/g and the coefficients Bq, Bi, • ■ • ,Cq^Ci^ ■ ■ ■ are obtained by the Fade 
coefficients in Eq.( |111.8]) . The degenerate solution of z is obtained when the line —sz becomes 
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tangent to the curve given by the r.h.s. of Eg. ( [III.9|) (Fig.2). This degenerate solution gives 
the critical point Zc- We apply various Fade methods for the r.h.s of Eq. (p.ll.9D . In Table II 
we list the critical point Zc as a Fade table form. From this table we see that the convergency 
of the Fade which involve even number of terms of original series is good. In the case that 
p < q, there appears a pole in the positive z and we consider the critical point is effected by 
this pole. We find the transition point Zc = 5.6, Sc = —0.4 and a\/g = lAhy [5, 3] Fade. 
In the one matrix model, we have an explicit expression which corresponds to Eq. (|III.9D 



as 



13 



-sz = —\U + 3z)2/3 + — + — (iii.io) 
108^^ ^ ^ 12 12 ^ ' 

and the exact transition point z = 4, s = —1/4 and a\/g = 16. The exponent 3/2 is of 
course related to the critical exponent of the free energy, 1 — 7st, where 7st is the exponent 
of the string susceptibility 7st = —1/2. In our case, this string susceptibility exponent is 



considered to be zero, since our system has a central charge c = 1 as shown in ref. [[13 
Thus, we expect that the singularity of the r.h.s. of Eq. (pi.9| ) at the negative Zc becomes 
(1 — z/zc) ln(l — z/zc). Then we must consider the logarithmic singularity and our simple 
Fade form lose the validity. However, we are interested in the region of positive Zc, and the 
assumption of Eq. ( [111.9]) may be valid. 

In the previous paper we discussed the flux lattice melting point from the direct 
calculation of the perturbation series for the GL free energy. In this study we used the 
convenient reduced temperature yt which is defined as 

yt= ■ (iii.ii) 



ieBl3/2'K 

This reduced temperature yt is expressed by the definition of z and s. 



Table III shows the obtained critical point in term of reduced temperature yt by the various 
Fade method. 

11 



For g < Qc, the saddle point Xc is frozen, and s becomes a constant Sc- The derivative of 
the free energy y = dF/dg is given by Eg. ( |111.7D as 



y 



2g 2g^ • 

Then the free energy becomes for the low temperature phase as 



;iII.13) 



Fjg) - FjO) 



SrCt 



- In Qf + 
2 ^ 2g 



(III.14) 



The phase transition occurs at 2; = = — Sea's /gc- 

In order to compare our model in the large limit with the usual = 1 case, we write 
the free energy F and the specific heat C in modified parameters in ref. Jl],^. Instead of 
g = j3eB / 4:Tia\ in Eq.( [11.6| ), we introduce a new variable g = j3eB / Ana^ , where a is a mass 
of Hartree-Fock approximation and related to as 



aB = a{l - Ag) 



(III.15) 
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peB 



(III.16) 



Then the Gibbs free energy G in two dimensions is given by 

G 



eB 
2^ 
eB_ 
2^ 



In 



a 



IT 



+ m 



(III.17) 



The function f{g) is equal to [F{g) — F{0)]/N'^ in Eg . (111.221 ). We have changed the variable 
g to g, since g is divergent at = 0. The new f{g) is readily obtained from Eq. (|II.22|) as 



r.2 , ".3 43^, 



-9 + -9 

2^ 9^ 



525^ 



6435^ 



133391 . 128589189285247 .0 
+1309 ^^_^_ / - 9471 / 



2263261- 



166966608033225- 



(III.18) 



This expression should be compared to the result for = 1 case; here we represent the 
previous result 
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2 29 

= -2~g -~g^ + A-f - 39— + 471.396594517^' 

Ar=l y oU 

:;;6 



6471.56257496^^ + ■ • • (III.19) 



It is easily seen that the N = oo case, the series is convergent, while the series for = 1 in 
Eq.( P^II.19| ) is an asymptotic expansion. 

For the low temperature region, we have from Eq.( |lII.l^ 



2 



m 



-hng + 8gsc (III.20) 



for g ^ oo. Since we have found Sc — —0.4, this behavior coincides with the low temperature 
free energy 0,^, of = 1 case, 

^ 'in^-i^ (III.21) 



( eB 
V 2tt 



2 " 1.16 

where a factor 1.16 is the Abrikosov ratio Pa =< > / < >^for the triangular 
lattice. 

The specific heat C, normalized by AC = v/(3 {v is volume), is obtained by the derivative 
of entropy S 

C 1 dS 



AC AC das 

and the entropy S is given by 



;ill.22) 



8C 

It may be interesting to compare the N = oo result of the specific heat with the previous 
result of = 1 case. The reduced temperature yt defined by Eq.( |11.17| ) becomes 

OB^ ^ (l-4g) 

\/¥ 

The specific heat in the high temperature region is analized by the Fade method, which has 
been studied before for = 1 case. Remarkably the large A^ result agrees with the previous 
result of A^ = 1 case. For the low temperature region, the specific heat is simply given by 
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C 1 



,2 



- 2sc (III.25) 



AC y2 

In Fig. 3, we represent the curves of the specific heat in the low temperature and in the 
high temperature regions. The values in the high temperature region are obtained by [5, 3] 
Fade result of Eq. (|III.9|) . In Fig. 3 , we represent also the curve of the specific heat of the 
= 1 case. There is a phase transition of third order at yt ~ —2.7 in the large limit. 

IV. DISCUSSION 

In this paper, we have developed a new series expansion for the matrix Ginzburg- Landau 
model and in the large iV limit, we have obtained a phase transition which corresponds to 
the superconductor transition in two dimensions. We have found the remarkable agreement 
about the specific heat between the large N limit and the usual GL model (A^ = 1) in the 
high temperature region and in the low temperature region. The phase transition point for 
N = oo case is higher than the value of = 1. We have obtained the transition point 
yt = —2.7, while we have yt = —10 for the usual Ginzburg-Landau model of A^ = 1. 

It is easy to evaluate the next order l/N"^, which corresponds with the diagrams of genus 
one. Also it is interesting to perform the numerical simulation by Langevin method or Monte 
Carlo method for the matrix Ginzburg- Landau model and to find the melting transition 
point which depends upon A^. We will represent this simulation result elsewhere. 

As a theoretical interest, our matrix model may be interesting in several points, (i) In 
one matrix model, the phase transition corresponds to the singularity of the density of state, 
where the density of the eigenvalue of the Hermitian matrix has a gap at the band center 



| 21| . Our gauged model is considered to be similar, although the eigenvalue representation 
is difficult due to the complex matrix. (ii)The matrix model in the large A^ limit represents 
the string behavior as a random surface in the double scaling limit at the negative g. We 
have discussed the phase transition at a positive g, and also the double scaling limit is 
expected at the phase transition point. In the superconductor, the vortex is a string, and it 
is interesting to think the phase transition, which we found in the large A^ limit, is related 
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to the string theory. We note that in a different context, the analogy of the phase transition 



in a strong magnetic field to string theory has been discussed ||2^ . As same as other matrix 



models [|T^,^, the transition may correspond to the condensation, like ideal Bose-Einstein 
gas. It is interesting to consider further the tachyon condensation for our phase transition. 
Our result of the new renormalized expansion becomes useful for the phase transition on a 
random surface in the case of the central charge c > 1. (iii) In the presence of the impurity, 
the melting phase transition is considered to be a second order |]25|j26[| , and a vortex glass 
or a gauge glass phase appears instead of the Abrikosov vortex lattice phase. This gauge 
glass phase has no long range order. It is interesting to note that for the matrix model in 
the large limit, below the phase transition point, there is no symmetry breaking, as seen 
in the density of state of the eigenvalue. The density of state has a gap, but still symmetric 
for the positive and negative eigenvalue. This is a common behavior seen in freezing phase 



transitions ||2^ . Thus a matrix model in the large A^ limit may become a model of a gauge 
glass state in the two dimensional superconductor in a magnetic field, in which the freezing 
transition is essential. 
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TABLE CAPTION 



TABLE I. Number of planar irreducible diagrams relevant in the derivation of Eq. (p.l3| ). 

TABLE II. Phase transition point obtained by various [p, g] Fade method in the < 
case. 

TABLE III. Reduced relative phase transition temperature which corresponds to the critical 
point Zc and obtained by [p, g] Fade method. 



FIGURE CAPTION 



Fig. 1 Flanar irreducible diagrams up to 0{g^). 

Fig. 2 Approximated [5,3] Fade form of r.h.s. of Eq. (|III.9|) is plotted against z by the 
solid line. The dotted tangent line determines the phase transition point Zc = 5.64 
and —Sc = 0.4 which is the slope of this line. 

Fig. 3 The specific heat scaled by AC against reduced temperature yt- For the derivation 
of this line f{g) is approximated by [5, 3] Fade form. Low temperature side is obtained 
from reexpressing Eq. (pi.l4| ) as the function of yt and evaluating its second derivative 
with respect to yt. The dotted line is the specific heat for the = 1 case obtained in 
ref. §. 
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TABLES 





Oifj) 


O(jf) 




0(jf) 




0(jf) 






# of planar diagrams 


1 


1 


1 


2 


3 


9 


22 


61 



TABLE I. 



q\p 


3 


4 


5 


2 


19.6 


5.64 


4.62 


3 


5.18 


4.12 


5.64 


4 


4.43 


5.19 





TABLE II. 



q\p 


3 


4 


5 


2 


-5.16 


-2.66 


-2.38 


3 


-2.54 


-2.23 


-2.66 


4 


-2.32 


-2.54 





TABLE III. 
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